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Abstract 

It is well established that at early times, long before the time of radiation-matter equality, the universe could have 
been well described by a spatially flat, radiation only model. In this letter we consider the whole primeval universe as 
a black body radiation (BBR) system in an n—dimensional Euclidean space. We propose that the (3 + 1) — d nature 
of the universe could be the result of a kind of thermodynamic selection principle stemming from the second law of 
thermodynamics. In regard the three spatial dimensions we suggest that they were chosen by means of the minimization 
of the Helmholtz free energy per hypervolume unit, while the time dimension, as it is well known was the result of the 
principle of increment of entropy for closed systems; that is, the so-called arrow of time. 


In 1989, Landsberg and de Vos [T] proposed a spatial 
n—dimensional generalization of the Planck distribution, 
the Wien displacement, and the Stefan-Boltzmann laws for 
black-body radiation (BBR) for a zero curvature space. 
Later, Menon and Agrawal [2] modified the n dimensional 
Stefan-Boltzmann constant found by Landsberg and de Vos 
by using the appropriate spin-degeneracy factor of the pho¬ 
ton without affecting normalized Planck spectrum given by 
Landsberg and de Vos. Shortly thereafter, Barrow and 
Hawthorne investigated the behavior of matter and radi¬ 
ation in thermal equilibrium in an n—dimensional space in 
the early universe, in particular they calculated the number 
of particles N, the pressure p and the energy density u [3]. 
More recently, Gonzalez-Ayala et al. [3] calculated several 
thermodynamic potentials for BBR such as the Helmholtz 
potential F, the enthalpy H, the Gibbs potential G and the 
entropy S by means of the generalized Planck distribution 
for an n—dimensional Euclidean space. Moreover, they cal¬ 
culated the corresponding densities per hypervolume unit 
for these potentials; that is /, h, g and s respectively. As it 
is expected, g and G resulted equal to zero for any spatial 
dimensionality n, such as it happens for n = 3. The rest 
of the potential densities s, /, and also the internal energy 
density u are functions depending on the absolute temper¬ 
ature T and the dimensionality n. When all these func¬ 
tions were plotted against T and n they seemingly did not 
show any interesting feature [3]. However, when a zoom was 
made in the region of very high temperatures (of the order 
of Planck’s temperature Tp) and low dimensionality (in the 
interval ne(l, 10)), a very surprising behavior was observed 
[3]. On the other hand, it is well established that at early 


times, long before the time of radiation-matter equality, the 
universe could have been well described by a spatially flat, 
radiation-only model [S]. However, in times very close to 
the Big Bang within the Planck scale, it is believed that the 
realm of string theory [6] or loop quantum gravity (LQG) [7] 
is found, where relativistic quantum gravity effects can be 
very important. In the present letter we analyze some ther¬ 
modynamic properties of a universe dominated only by radi¬ 
ation within a flat spatial geometry, beginning at the Planck 
scale and until times around t = lO^^s, the time separating 
the radiation-dominated epoch from the matter-dominated 
epoch [8]. Our analysis is based on an n-dimensional Eu¬ 
clidean space filled with BBR. As it was shown in [3], this 
kind of system exhibits convexities for certain thermody¬ 
namic potential densities only in the region of very high 
temperatures and low dimensionality. In particular, the in¬ 
ternal energy density u, and the entropy density s show in 
the first place (as n increases) a local maximum and then 
a local minimum while in the Helmholtz free energy den¬ 
sity first a local minimum and then a local maximum. It is 
very interesting that the first function that presents a criti¬ 
cal point is the Helmholtz potential density. These extreme 
points are found within the context of isothermal processes 
in a plane of the corresponding potential density against the 
dimensionality n. Remarkably, the form of the isotherms is 
reminiscent of the isotherms found in some first order phase 
transitions in gas state equations of the Van der Waals type, 
for example (see Fig. 2). If we start at Planck temperature 
Tp, the density / finds its minimum value at T* = 0.93Tp, 
and the other two densities {u and s) present their extreme 
points after this value. If n is considered as an integer, the 
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minimum of / is located at u = 3. For any later moment, 
after the time corresponding to T* = O.QSTp, the spatial di¬ 
mension n = 3 remains “frozen”. This extreme point occurs 
in the unique isotherm corresponding to a saddle point in 
the plane / vs n. Below this point a continuous transition 
from n = 3 to any n 7 ^ 3 is forbidden by the second law of 
thermodynamics and for all T > T* = 0.93Tp, any value of 
n is permitted (allowing for proposals of the type of string 
theories [ 0 ] or even of the type of the so-called vanishing di¬ 
mensions models [nmniiii]). This kind of thermodynamic 
behavior offers a possible mechanism for determining the 
3-1-1 dimensional nature of the space-time. 

The question of why is space 3—dimensional goes back to 
ancient Greece |12j . In modern times this question was first 
raised by Kant in 1746 m- Later, Ehrenfest in 1917 by 
means of the stable orbits argument showed that n = 3 
M- In 1983, Barrow brought forward a very interesting 
approach to the dimensionality problem |12| . Since then, 
many authors have worked on this problem extended to the 
case (3-1-1); that is, including time. Such is the case of 
Brandenberger and Vafa ng, that in 1989 proposed a nat¬ 
ural mechanism for explaining why there are 3 large space 
dimensions in the context of string gas cosmology ngns]. 
Within the context of the (3-1-1) problem Tegmark pub¬ 
lished an enlightening article summarized through his Fig. 
1 HZ]. For a deeper discussion on this issue one can also see 
the following works [TH][TT] [ini 1^ . 

For a black body in an n-dimensional space it is known 
[H El EHH EH [ 22 ] that the number of modes per unit fre¬ 
quency interval is, 

2 (n-l) 7 rtF:/”-i 
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where V is the hypervolume and n the dimensionality of 
the space. The multiplication of the above equation by the 
Bose-Einstein factor, which is the probability of finding a 
photon with a frequency v gives the number dN of photons 
in the frequency interval v to v + dv, 
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respectively. Since the Helmholtz potential density will be¬ 
come very important in our analysis, the complete form will 
be given [4], 

, 2(n-l)^t(A:T)"+iC(n + l)F(n + l) 
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The corresponding thermodynamic potential densities are 
immediately obtained by dividing each quantity of Eq. m 
by the hypervolume V. All these results agree with the well- 
known 3-dimensional cases |23] . 

In a previous work |4] we showed that the thermodynamic 
potential densities have regions of minima and maxima when 
we consider isothermal processes. In Figure 1 the critical 
points of u(n^T), f (n,T) and s(n,T) with respect to the 
dimension n {dnU, /, s = 0) are shown. These critical points 
can be minima, maxima or saddle points. Notice that for 
each function there is a certain critical temperature T* at 
which for any T > T* there are no more extreme points 
of the corresponding thermodynamic potential density. The 
region of transition between such convexities (from maxima 
to minima or viceversa) occurs at dimensionalities in the 
interval from 3 to 5 and at very high temperatures (nearly 
the Planck temperature). In each case (see Fig. 1) the 
boundary between the maxima and the minima is the sad¬ 
dle point corresponding to the unique temperature T*. In 
Fig. I it can be seen that by decreasing the temperature 
from T = Tp, the first function that reaches a critical point 
is / (n, T). The isotherm that passes by this point is that 
with T* « 0.93Tp. For temperatures below this value there 
are maxima and minima of the Helmholtz density function 
and it is possible to talk about thermodynamic optimization 
criteria. That being said, the first function to optimize was 
the Helmholtz free energy density and it happens in a region 
near the value n = 3. This might not be a mere coincidence 
as will be pointed out later. 
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The integration of hvdN/V for all the frequencies gives the 
energy density (per hypervolume unit) u(T,n), which is 
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Figure 1: -^f,s,u = 0. The regions of maxima and minima 
are well defined. The red continuous lines are the minima 
and the dotted lines are the maxima of /, s and u. The blue 
vertical dashed lines are located at the unique saddle point 
of each potential. 


The Helmholtz function F, the entropy S, the pressure p 
and the Gibbs free energy G are given by the following ex¬ 
pressions |3] 


F = -—u, S = 
n 



Vu 


T 


P = 


(1)„,G = 0. 


(3) 


The separations between the regions of maxima and minima 
for u (n, T), / (n, T) and s (n, T) are well dehned (see Fig. 1 
and it is also summarized in Table 1 of Ref. HI)- 
Since the works by Kaluza and Klein |241 E5] . many pro¬ 
posals about universe models with dimensionality different 
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from three have been published [U [21 [3J [HI UHl (HI IHl IHl 
1211 mi 1211 mi I2H1 mi uni mi • Remarkably this has been 
the case of results stemming from the string, D-branes and 
gauge theories dSl lISlEelETlEllEalSQ]. However, nowa¬ 
days we only have evidence for a universe with three space 
and one time dimensions. 


Temperatures as high as are only possible in very 

early times in the evolution of the universe. It is known 
that this period was dominated by energy in the form of ra¬ 
diation, which means, that in this period the universe could 
have been a flat space [32] • Thus, considering the whole 
primeval universe as a black body radiation system in an 
Euclidean space is in principle a reasonable approach. It 
is considered [331 153] that the critical energy density at the 
Planck epoch (t Itp ss a Planck time) was around 

u Ri lEp/Vp ft! 10® J/Vp ( an energy of Planck in a Planck 
volume) and the temperature around T « ITp Ri lO^^K 
. The temperature and energy density obtained when the 
maximum of the Helmholtz density is located at n = 3 (in 
fact n Ri 3.3, see Fig 2) are surprisingly close to these val¬ 
ues. For example, according to string theories at some point 
(at the end of the Planck epoch) the rest of the dimensions 
collapsed (or they simply stayed at the same size) and only 
the 3—dimensional space grew bigger |6l|32]. The remaining 
question is: why did only 3 dimensions expand? If thermo¬ 
dynamic laws were born with the universe, then, thermo¬ 
dynamic analysis could give us a clue that a maximum or 
a minimum criteria would have been fulfilled at the earliest 
period of our universe. The conditions of energy and tem¬ 
perature in that epoch were probably suitable to maximize 
some kind of thermodynamic function. 


In Figure 2 we depict the isotherm curves corresponding to 
the negative of the Helmholtz potential density. In fact, from 
Figure 1, as the temperature decreases from T = ITp to 
zero, the first potential that reaches a critical point (a saddle 
point) is the Helmholtz energy density. The first isotherm 
that reaches a critical point is the one with T* = 0.93Tp. 
Between T = ITp and T* there are no other critical points. 
After this moment, that is, for temperatures in the interval 
from T < T* to T = 0, there is always first a maximum 
and then a minimum of this function. This kind of behav¬ 
ior resembles the form of the isotherms of a gas-liquid phase 
transition. As is well known the critical isotherm (Cl) corre¬ 
sponds to a saddle point which divides the pressure-volume 
plane in two regions: above Cl, the transitions do not occur 
and below Cl, the transitions are permitted. However, as we 
shall see in the present case, in the region below the criti¬ 
cal isotherm T* restrictions imposed by 2nd law appear and 
transitions between dimensions are forbidden. With this in 
mind, in Figure 2 we can see the isotherms as possible tran¬ 
sitions in the space dimensionality. 



Figura 2: The continuous line are the isothermal processes, 
and in dotted line the critical points. The highest extreme 
point corresponds to the saddle point located at the isotherm 
T = T*, above it there are no more extreme points. 


It is very interesting to observe that the critical isotherm T* 
in Fig. 2 divides the plane f — n in two regions; above T* 
there are no restrictions over any particular value of n. How¬ 
ever, as mentioned before, below T* there will be restrictions 
on transitions from n = 3 to any other dimensionality. 


A well known theorem derived from 2nd law linked with the 
Helmholtz potential F is the Helmholtz potential minimum 
principle [35] . This theorem implies that for isothermal pro¬ 
cesses 




T 


< 0. 


(4) 


where the equality is fulfilled by reversible processes. Let’s 
consider a hypercube with volume V = i?" as the black 
body system, where R is the length of the edge. The above 
equation can be expressed as. 
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dfV 

dV 


^ dndV - ^ 


where it has been used the consideration that §p = 0. Then, 
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In an adiabatic expansion (no heat exchange between the 
universe and the “exterior”, a reasonable assumption already 
considered in standard cosmology), the condition showed in 
m is = constant. By using Eq. Q and the ex¬ 

pression for p in Eq. the adiabatic condition for two 
different moments is 


R = 


TqRq 

T 


Tplp 
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( 6 ) 


where Rq and Tq are some initial conditions that, in our 

case, are very near the Planck epoch and Ip Ri 1.6 x 

is the Planck length. This is a result that agrees with typical 
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treatments in cosmology (see for example [33]). Then, Eq. 
((Sj) is now (in Planck units) 

( 7 ) 


persion relation given in |7|, obtaining the following result 

01 , 

,, 2^t(n-l)T"+i fr(n + l)C(n + l) , 

/ (n, T) = --+ 

r(f) ( n 


The difference between the dot-dashed line in Fig. 2 and 
the restriction given by Eq. ([7|) (red dotted line in Fig. 2) is 
very small. According to Eq. ([7|) for an isothermal process 
the region below the dotted line is prohibited. That is, once 
the integer value n = 3 is reached at T*, for any T < T* 
the space dimension of the BBR system remains frozen at 
this dimension because, as it can be seen in Fig. (2), any 
transition between n = 3 and n 3 is not permitted by 2nd 
law. 

As mentioned before, the first objective function reaching an 
extreme point was / (n, T), but is it an adequate quantity to 
optimize? Notice that when the internal energy is minimized 
it is possible to have an infinite number of configurations, 
some with more order than others for the same value of the 
energy. On the other hand, by maximizing the entropy it 
is possible to obtain an endless number of energetic config¬ 
urations, thus, by themselves neither the maxima nor the 
minima of s and u, respectively, determine an advantageous 
and unambiguous objective function. This is not the case 
for f = u — Ts, which gives a kind of trade off between en¬ 
tropy (organization) and energetic content, and this could 
be indeed a more meaningful optimization criterion. When 
maximized, it offers the better configuration that gives the 
less amount of disorder for a given energy. Lets recall that 
in stability analysis for closed systems, the stability points 
are found in the maximum value of entropy or in minimum 
energy configurations. Each one on their own do not es¬ 
tablishes an optimum criterion, because a restriction over u 
(or s) allows an infinite number of configurations of s (or 
u, respectively). Then, the minimum of the Helmholtz free 
energy gives a good commitment between the organization 
of a system and its energetic content. In the scenario where 
the universe cools down through an adiabatic free expan¬ 
sion, the first potential whose optimization is reached is the 
—/ > 0 function. The first maximum appears at u = 3. 
We emphasize the fact that the only isotherm with a saddle 
point {dnf = 0 and dnnf = 0) is that with T = .9STp, 
which could have fixed globally or locally the dimensional¬ 
ity of the black body system to n = 3. After this point, the 
adiabatic cooling process continues due to the expansion. 
Nevertheless, any isothermal dimensionality “phase” tran¬ 
sition is forbidden as soon as the temperature diminishes, 
because for small times St after, the change from n = 3 to 
any other value of n is inside the forbidden zone. It is possi¬ 
ble to see in Fig. l.b that the saddle point of / (n, T) is not 
located exactly in n = 3. This could be a consequence of 
having a simplified model. A possible correction using the 
results of quantum gravity can be incorporated by means 
of the modified dispersion relation |7] stemming from LQG 
theories. In ref. |4] we found the modified Helmholtz free 
energy density / by taking into account the modified dis¬ 
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where a and a' take different values depending on the de¬ 
tails of the quantum gravity candidates. By manipulating 
the value of the parameters a and ol in Eq. ([5|), which de¬ 
pend on the LQG model, it is possible to modify the value at 
which the saddle point occurs. In Fig. 3, it is depicted the 
analysis of maxima and minima for the Helmholtz potential 
density incorporating the LQG corrections. A similar anal¬ 
ysis of critical points reveal that the corresponding saddle 
point T* decreases in one order of magnitude. 



Figure 3: Maxima (dotted lines) and minima (continuous 
lines) of / incorporating LQG corrections. In the case 
a = 0.1 and a' = 2.28 (these values reproduce the curves of 
the spectral energy density showed in H) the saddle point 
is located at n ~ 2.44. In the case a = 0.15 and a' = .001 
the saddle point is at u 3 (see vertical shaded lines). 


In this way it has been proposed a model based on simple 
suppositions to study scenarios that put the space dimen¬ 
sionality n « 3 as a convenient candidate to optimize a 
thermodynamic quantity: the Helmholtz potential density. 
After picking this value, later transitions were prohibited 
by second law restrictions and the dimensionality remained 
“frozen” at n = 3. In regard the 1 in (3 -I- 1) — d, is generally 
accepted that the arrow of time was imposed by the princi¬ 
ple of entropy increment for closed systems. Thus, possibly 
the laws of thermodynamics selected the dimensions of the 
universe. 
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